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Abstract 

A thermodynamically consistent framework able to model either diffusive and 
displacive phase transitions is proposed. The first law of thermodynamics, 
the balance of linear momentum equation and the Cahn-Hilliard equation for 
solute mass conservation are the governing equations of the model, which is 
complemented by a suitable choice of the Helmholtz free energy and consis- 
tent boundary and initial conditions. Some simple test cases arc presented 
that demonstrate the potential of the model to describe both diffusive and 
displacive phase transitions as well as the related changes in temperature. 
Keywords: diffusive and displacive phase transition, thermodynamics, 
phase field model. 

1. Introduction 

Solid-to-solid phase transitions can be of different nature: diffusive and 
displacive. Diffusive transformations are second order phase transitions, and 
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are able to describe, among others, processes involving diffusion of atoms. 
Typical examples include spinodal decomposition, vapour-phase deposition, 
crystal growth during solidification and grain growth in single-phase and 
two-phase systems. Spinodal decomposition can occur, for instance, in an 
eutectoid steel: a diffusive transition from a high temperature stable phase, 
austenite (Figjl^), to perlite (Fig|l|D) takes place at low cooling rate; as the 
temperature decreases slowly, carbon (the solute) can migrate outside iron 
cell and originate cementite layers. 

Displacive transformations are first order phase transitions involving a 
lattice distortion. Typical examples include the martensitic transformation of 
shape memory alloys or martensitic transformation of quenched steels: here 
the displacive transformation from austenite to martensite (Figjlj:;) occurs at 
high cooling rate; the temperature decreases fast, and there is no time for 
diffusive phenomena to take place so that the new phase originates with a 
lattice deformation. 

Diffusive and displacive phase transitions are concurrently involved in 
many and very important industrial processes, such as the heat treatments 
of steel. In particular, in a steel subjected to a heat treatment there may be 
different zones in which, due to a different cooling rate, the transformations 
are of different kind, leading to different phases. 

A number of approaches may be adopted to model both the aforemen- 
tioned phase transitions; for instance, use of phenomenological laws can 
be made: Brokate [2J, following the approach proposed by Hoemberg [3], 
presents a model constituted by the Koistinen-Marburger rule for the pre- 
diction of the martensitic phase fraction, the Johnson-Mehl formula for the 
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estimate of the perlitic phase fraction, together with the Scheil's additivity 
rule for non-isothermal conditions. However, in Author's opinion, such an 
approach does not make clear the physics behind the phenomena involved. 
A different way to tackle the problem may consist in setting a framework 
in the context of nonlinear thermoelasticity including the Ginzburg-Landau 
theory of phase transitions [5]. Following the phase- field approach, a 
microstructure is identified by a variable named order parameter. The phase 
transitions rely on phase evolution equations; constitutive laws are expressed 
as the derivatives of the Helmholtz free energy density, which is given as a 
function of the order parameters. 

In particular, a Ginzburg-Landau theory was proposed by Falk [B] for the 
martensitic phase transition in shape memory alloys using an order parameter 
dependent on the strain tensor. The same idea has been developed by Barsch 
et. al. [7] to describe the microstructure of inhomogeneous materials and 
their proposed model has been extended to simulate martensitic structures 
in two- and three-dimensional domains. In the work by Ahluwalia et. al. [8] 
this theory is generalised to describe a two-dimensional square to rectangle 
martensitic transition in a polycrystal with more than one lattice orientation. 
Other models which describe the martensitic transformation in a Ginzburg- 
Landau framework using a different order parameter not dependent on the 
strain tensor are proposed by Levitas et. al. 0, Wang et. al. [lU], Artemev 
et. aim], Berti et. al.[I2]. 

Ginzburg-Landau frameworks have also been applied to study diffusive 
phase transitions. In this context, Cahn |T3] presented a first model in the 
simplest case of diffusive phase transition without accounting for thermal 
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effects. An extension to Calm's model is presented - among otliers - by Alt 
et. al. [H]; in their paper the coupled phenomena of mass diffusion and heat 
conduction in a binary system subjected to thermal activation have been 
modelled. Other models accounting for the mechanical aspects related to 
the diffusive phase transitions are the ones by Onuki et. al. [15] and Fried 
et. al.dS]. 

First attempts to study diffusive and displacive phase transitions together 
were made by Rao et. al. and by Levitas [IB]; the same idea has 

been developed by Bouville et. al. [19], who proposed a Ginzburg-Landau 
framework including the balance of linear momentum equation and the Cahn- 
Hilliard equation and analysed the effects of a volume change consequent to 
a diffusive or a displacive phase transformation. 

In this paper, an attempt to describe both displacive and diffusive phase 
transitions in a thermodynamically consistent framework is made; the cou- 
pling between thermal, chemical and mechanical effects is accounted for, and 
a way to overcome the difficulties arising from the treatment of the gra- 
dient terms is proposed. Three are the governing equations, presented as 
balance equations for observable variables: the balance of linear momentum, 
the Cahn-Hilliard equation as a solute mass balance [25] and the heat equa- 
tion (balance of internal energy). The model is completed with a suitable 
description of the free energy. The analysis is reduced to a two dimensional 
setting, which is simpler than a three dimensional one, but still meaningful. 

The outline of the paper is as follows: after explaining the choice of the 
order parameters for each transformation in Sect|2| in Sect j3] balance equa- 
tions are provided for linear momentum, solute mass and internal energy; the 
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(a) (b) (c) 

Figure 1: Microscope image of: (a) austenite, (b) perlite and (c) martensite [T]. 



first two will act as the evolution equations needed for a Ginzburg-Landau 
framework. In Sect|4]the restrictions provided by the second law of thermo- 
dynamics are highlighted and its exploitation in the form of the Clausius- 
Duhem inequality is provided. In SectjS] the model is completed with the 
description of the Helmholtz free energy density, which will allow an explicit 
form for the constitutive laws to be derived in SectjSj The resulting model 
is a diffuse interface one, and this implies the presence of the gradient of the 
order parameters, which will put some issues regarding the exploitation of 
the second principle of thermodynamics (Sectjij) and the necessary high or- 
der boundary conditions, as highlighted in Sect |6.3[ The paper ends drawing 
some conclusions. 

2. Model of the problem 

The model is set at a microscopic scale of observation which may be, for 
instance, a single-grain scale and relies on a number of order parameters to 
describe the different kinds of transformations involved. In order to account 
for the presence of three phases (austenite, perlite and martensite), two order 
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parameters playing as the unknowns of the problem may be introduced. 



2.1. Diffusive order parameter 

To account for diffusive phenomenon, an order parameter called c is set: 



being Vt a two-dimensional domain. 

This parameter describes the evolution from one phase to another one 
with a different structure. According to experimental observations and the 



represent the solute mass fraction. Hence, Pearlite, and in general phases 
obtained by diffusive transformations, will be the one for which this order 
parameter is non-zero, i.e. the carbon has migrate to generate areas in which 
there is a great quantity (c > 0) in spite of other areas, in which there is a 
lack of carbon (c < 0). Moreover, the higher value c will assume in a point of 
the domain, the more the transformation will be at a late stage there. Given 
this interpretation, the Cahn-Hilliard equation may be seen as a balance of 
mass equation. 

2.2. Displacive order parameter 

Regarding martensitic phase formation, a shear strain needs to be applied 
to the lattice in order to obtain the phase change [20], [21]. For this reason, 
in considering displacive phenomena, the common choice for the order pa- 
rameter, called 62, in a two-dimensional setting [S] [22] [23] is 




(1) 



use of the Cahn-Hilliard equation |25j (see Sect|3|, the order parameter may 



£^11 — ^22 



(2) 



62 = 



V2 
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where 

e = 1 (Vu + Vu^) (3) 
is the hnearised strain tensor and 

w = w(cc,i), for(a;,t)ef2 x [0,oo) (4) 

is the displacement field. 

3. Balance equations 

We start the description of the model by stating the balance equations 
and then by examining the restrictions imposed by the second law of thermo- 
dynamics. We will last enter appropriate constitutive hypothesis by choosing 
a constitutive law for the free energy. 

The model consists of three balance equations, reported here in their 
pointwise formulation: 

• balance of linear momentum: 

pii - V • o- = p/, (5) 

where a is the stress tensor and / is the external body force; 

• balance of solute mass: 

c = -V • J, (6) 

where J = cV is the flux vector related to the solute mass fraction {V 
is its the velocity) ; 
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• balance of energy: 

e = pin + pl + h; (7) 

here, e is the internal energy density, pl^ is the mechanical internal 
power (density), pi is the internal power associated to the solute mass 
balance equation and h is the internal thermal power. This can be 
expressed by the thermal power balance 

/i = - V • q + r, (8) 

q being the heat flux and r the external heat supply 

The balance of linear momentum and of solute mass give rise to two power 
balance equations. 

To account for non-local effects and spatial variations of the order pa- 
rameters (interfaces), we demand a non-local constitutive relation for stress 
and chemical potential. Some cares are needed due to the presence of non- 
localities in the constitutive equations, as this will lead to non-conventional 
expressions for the internal powers associated to momentum and solute mass 
balance. 

For the time being it suffices to consider this relation in the quite generic 
form 

<T = <T;oe(r;oe,r;oc)-v-i](r), (9) 

where Tioc — (e, c, T) are the local state variables and F = {^loc: Ve, Vc) 
is the full set of state variables Moreover, (Tjoc, V • S are supposed to be 
symmetric second order tensors. The dependence upon Tioc accounts for 
possible dissipative contributions to the local stress; we will not consider 
dissipative contributions to the non local part of the stress. The balance of 

8 



mechanical power is obtained by multiplying the balance of linear momentum 
equation hj v = ii, which leads to 

iQp^') =[V-(^.oc-V-I]) + /].« = 

= -{(Tioc ■■ i + SlVe) + V ■ [cr • V + S : s] + / ■ t;, (10) 
having enforced the symmetry condition for the tensors and therefore sub- 



stituted £ with Vv. We interpret (10) as a balance equation for the rate of 



change of kinetic energy and the internal and external (mechanical) powers: 

EiS) = I ^-pv'dw, (11) 

KXS) = j{<Tioc ■■ e + s;Vs) dw, (12) 

^m('5) = / f-vdw+ [ (cr • V + S : 6) ■ i^ds, (13) 
Js J as 

where S is any sub-body, dw is the volume measure, ds the surface measure 
and u the surface normal. We have used capital letters for the integrated 
powers, so that Vl^{S) = f^pl^dw. 

The balance of mechanical power equation is expressed as 

E{S)+VUS)=V^^{S). (14) 

Now we consider the power associated to the balance of solute mass equation. 
We call /i the chemical potential, which is the variable conjugated to c defining 
the power balance associated to 

Hc = -/iV ■ J. (15) 

A constitutive expression function of F for the chemical potential has to be 
given, as well as a non-local relation: 

fi = fiioc{Tioc)-V -MiT). (16) 
9 



We can rewrite (15) as 

I2{C + V ■ J) = CIJ. + V ■ (fiJ) -J-Vi2 = 

{c^Mioc + M -Wc- J ■ V/i) + V ■ (/i J - cM) = 0. (17) 

We define the internal and the external powers as follows: 



Vl{S) = / (cM - fiJ) ■ uds 
Jas 

so that the balance of powers in the sub- volume S is expressed as 

Vl{S)=V::{S). (18) 

Note that there is no energy related to a second order time derivative in the 
balance of solute mass equation. 

4. Restrictions from the Second Law of Thermodynamics 

We will examine the restrictions entailed by the second law of thermody- 
namics in the form of the classical Clausius-Duhem inequality: 

77>-V-| + ^. (19) 

By using the first law of thermodynamics (balance of energy equation): 

e = - V ■ g + r (20) 

(where = pl^ + ) to eliminate the source term r, and introducing the 
Helmholtz free energy density = e — Trj, the Clausius-Duhem inequality 
can be expressed in the more convenient form {reduced inequality) 

i^ + f^-vlot<-^^T. (21) 
10 



By substituting the expressions for the internal powers and remembering that 
■0 is a function of the variables F = (e, c, T, Ve, Vc), the inequality becomes 

(Vt + v)T + {ipc - l^ioc)c + (V'vc - M) ■ Vc+ 
+ (^ve - S)!V£ + {i/Je - CTioc) : £ + J ■ V/i + I • VT < 0, 

where the subscripts represent partial derivatives. This inequality has to be 
satisfied for all processes (e, c, T, Ve, Vc, V/x, VT). Since r] (as well as ■0) 
depends only on the state, we obtain 

V = -^T, (22) 

by considering processes with only T different from zero. Similarly, as also 
/iioc, M-i S are assumed state-dependent only, we have 

tJ'ioc = V'c, (23) 
M = ^vc, (24) 
S = V'v.. (25) 

The entropy inequality is therefore reduced to 

[V'e-o-^oc] :e + J-V/x+|-Vr<0. (26) 

The current J and the heat flux q are in general functions of the state and the 
process. For our purposes it will be sufficient to make the following simple 
assumptions ensuring the exploitation of the inequality: 

CTloc = V'e + fan, CTan^l £, (27) 

J = -aV/x, a(F) > 0, (28) 
q = -kVT, k{T) > 0, (29) 
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where 7(r) is a positive definite fourth order tensor of which we provide an 
exphcit formulation in Sect. |6.1[ a and k are positive coefficients in general 
dependent on F, more often only on the temperature T. The subscript in 
the symbol cr^^ stands for anelastic. The presence of a dissipative term in 
the local stress provides a damping mechanism for the displacive transition, 
as showed in Sect. ??. 

We point out that the chemical potential and the non-dissipative part of 
the stress can be written accounting for their local and non-local parts in 
terms of functional derivatives of a free energy functional 

m= [ ^{T)dw, (30) 



where VL is the material domain and dw the material volume measure. Then, 

/i = V^c - V ■ V'Vc = Cr - (Tan = z/'e - V ■ z/'Ve = (31) 

OC OS 

these definitions being consistent with the choice of the boundary conditions 
made in Sect J6.3[ 

It is possible to substitute the constitutive relations in the equation for 
the internal powers, obtaining the following internal power density: 

ftot = {^e + (Tan) : i + z^vs^Vs + + " Vc + a|V/ip = (32) 

= ^-TV'T + a|V/i|' + £: (7:s). (33) 

Besides, by using r] = —ipT, we have 

e=j^{ilj + Tr,)=^-t^T- T^fr, (34) 



so that (20) can be rewritten in term of the free energy as 

- TiPt = a|V/i|^ + £ : (7 : e) + V ■ (kVT) + r. (35) 
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This is a thermodynamically consistent framework which will be substanti- 
ated with a convenient choice of the free energy %jj (Sectlsl). 



5. Free energy 

Since constitutive equations have been determined by means of derivatives 
of the free energy, the choice of this latter is a convenient way to characterise 
the constitutive content of the model. There are two main phenomena in- 
volved in the discussion: the diffusive and the displacive transitions, which 
have to be accounted for at the same time; this is done by exploiting the 
additivity property of the energy. 

The free energy density is written as the sum of different contributions: 

ip = i)o{T) + i)thei{c, 62, T) + i^diffic, T) + i^dispiie, T) 

+ 1pcpl{c, 62) + V3rad(Vc, Vca), (36) 

where: 

• V'o is a term accounting for the background thermal properties of the 
material (contributing to specific heat); 

• 'ipthei is a term accounting for the dilatation due to phase transitions 
and for the contribution of the shear strain; 

• i^disph which is a function of 62, accounts for the displacive transition; 

• i^diff, which is a function of c, accounts for the diffusive transition; 

• ipcpi is a coupling term for the interaction between the two phases; 
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• "^grad IS a gradient term and accounts for non-local properties and spa- 
tial variations of the order parameter. 

In the Ginzburg-Landau theory of phase transitions, it is assumed that the 
free energy density admits a power series expansion with respect to an order 
parameter. 

5.1. Diffusive part of the free energy 

The diffusive part of the free energy may be postulated as follows [T3] : 

4 2 

i^d^ff = A^ + BiTf-, (37) 

where A is a constant. The function B{T) is given by: 

BiT) = 5o^^, (38) 
J-p 

with Bq constant and Tp the temperature of perlitic phase transition. In 
Figj2| ipdiff as a function of c is depicted for various values of temperature. 
At high temperature there is only one stable phase identified as the minimum 
for c = 0, while at low temperature a different stable phase appears for values 
of c starting from for T = Tp and increasing for decreasing temperature. 
The transition occurs without hysteresis, i.e. with a continuous variation of 
the order parameter. 

5.2. Displacive part of the free energy 

The displacive part of the free energy may be postulated as follows [6], 
123: 

^Pd^s,l = D'|-E'j + F{Tfj, (39) 
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Figure 2: Diffusive free energy. 

being D, E constants. The function F{T) is given by: 

T-Tm 



F{T) = Fo 



T. 



M 



(40) 



with Fo constant and Tm the temperature of martensitic phase transition. 

FigjS] shows that at high temperature only one phase is stable, identified 
with the minimum for 62 = 0. For T = Tm another phase become stable, 
as the free energy shows two minima. Below the critical temperature, the 
phase corresponding to the order parameter equal to zero is unstable, whilst 
the other one, corresponding to non-zero values of 62, is stable. Note that, 
as 62 describes a first order phase change, the transition between the phases 
occurs with a jump in the order parameter. 
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Figure 3: Displacive free energy. 

5.3. Thermo- elastic part of the free energy 

The expression of ipthei is the one for the free energy accounting for the 
dilatation due to phase transitions [B] [21], changes in temperature and for 
the contribution of the shear strain. 



being ei 



1 



(ei - xicC - xucl - (5(T - Tq))' 



-3; 



72 



^11 + ^22), 63 = £12, G, H, To, Xic and X12 constants. 



(41) 



5.4. Coupling part of the free energy 

ipcpi is the couphng free energy and penahses the formation of one low 
temperature stable phase (e2 7^ or c 7^ 0) with respect to the other [6] [21]. 



(42) 



being X2c a constant. 
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5. 5. Gradient part of the free energy 

The gradient part of the free energy is chosen to be the simplest possible 
one [6] [2lj: 

Ipgrad = ^e^^ + ^c^, (43) 

with Kf. and constants; it accounts for non-local effects and spatial vari- 
ations of the order parameters (interfaces) . 

6. Summary of the equations of the model 

6.1. Constitutive relations: stress and chemical potential 

The constitutive laws of stress and chemical potential can be put in an 
explicit form. It is useful to introduce the following notation: 

ei = ^(*®* + j® j), (44) 

^2 = 3 ® j). (45) 

ez = i®j+j®i. (46) 
The stress tensor can now be written as 

cr = (Tio, - (KeV'e2)e2, (47) 
where the local part is given by four contributions: 

C/oc = displ + (^thel + cpl + (^ani (48) 

with 

CTdispi = (^De^' - ^62=^ + Fo ^'J^' eaj e^, (49) 

(Tthei = G (ei - xicC - xi2el - /3(T - To)) (ei - 2x126262) + He^e^,, (50) 
cTcpz = 2x2cC^e2e2, (51) 
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as it regards the conservative contributions; for the dissipative term, we 
assume 

cTan = 7e2e2, (52) 

which amounts to consider 7 = 762® €2. In this way, we consider a dissipative 
term associated only to the displacive order parameter 62; its effect is to dim 
the time oscillations in the order parameter after the displacive transition. 
The chemical potential is the following: 

= fJ'ioc - KcV'^c, (53) 

where 

filoc = Ac^ + Bq—:, C 

-Lp 

- Gxic{ei - xicC - Xi2e2^ - f3(T - Tq)) + 2x2ce2^c. (54) 
These expressions have to be substituted in the first two equations of system 



(55). 



6.2. Equations of the model 

We now summarise the equations of the model. We collect the three 
balance equations: 

c = V ■ (aV/i) 

pii = V ■ (7 + b 

w \ / \ 1 (55) 

c,(T)T - T yi^^c + G/3a;iej c + (^^ + 2G/3xi2 j 6262 - G^e^ 

-a| V/ip - = V ■ (kVT) + r, 
where the last one is the complete heat equation, obtained working out the 



expression of ifjx which appears in eq. (35) and accounting for the anelastic 
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stress (Tan in (52). We remark that 

c,(T) = -Ti^TT = -T{ {^o)tt + G(3^ ) (56) 

denotes the specific heat. The given expression of the heat equation allows 
to identify the thermo-chemical and the thermo-mechanical couphng terms 



(contained in the second term of eq. 55}), which contribute to the changes 



in temperature during a phase transition, and the chemical and mechanical 
dissipations, respectively a|V/xp and 763, which act as heat sources when a 
displacive or a diffusive transition occurs. 



Note that the first of ( 55 ) is the classical Cahn-Hilliard equation 



6. 3. Boundary and initial conditions for the differential system 

Well-posed boundary and initial conditions are necessary to complement 



Eqs.(55). Let f2 be the domain in the material description (for the purposes 
exposed here, no distinction is made between the material particles of the 
material description and the reference configuration of the referential descrip- 
tion). To state the boundary conditions in a general mixed form, we assume 
that the boundary dfl is partitioned into two subsets in two different ways: 

dn = Oils U d^u = d^g U d^T, (57) 

with 

dn, n driu = 0, dng n = 0. (58) 
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We assume the following boundary conditions: 



Vc(a3, t) ■ V 


dn 


-- 


V/i(a;, t) ■ V 


dn 


= 


Ve2( 


X, 


t) ■ V 


an 


= 


q ■ V 






0-loc{x, 


t) ■ V 




= t 








-u{ 


X) 


T{x, 


t) 


dflT 


= f 


[x). 



(59) 



tix) 



We remark that the first and the second conditions ensures the global con- 
servation of solute mass in the domain: 



d 
dt 



cdw 



J ■ ly ds = 0. 



(60) 



an 



The non-local boundary conditions (Eqs.(59)2, (pOjls) are needed in this 
model due to the non-local form of the stress [26]. The other conditions 
are classical. 

The initial conditions are the following: 



u{x, 0) = Uq{x) 

u{x, 0) = Vo{x) 
c{x,0) = Co{x) 
T{x,0)=To{x) 



(61) 



7. Conclusions 



A model to study both diffusive and displacive phase transitions has been 
proposed in a thermodynamically consistent framework in which, for the 
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presence of the balance of energy equation, non-isothermal conditions can 
be accounted for. The resulting model is a diffuse interface one, and re- 
lies on the choice of two order parameters and of its gradients; this rises 
some issues about the exploitation of the Clausius-Duhem inequality, and 
a way to deal with this has been proposed. In particular, we have consid- 
ered a suitable modification of the first law of thermodynamics by including 
non-conventional expressions for the internal powers (due to non-localities), 
keeping at the same time the classical Clausius-Duhem form for the second 
law of thermodynamics. 

The equations of the model are the balance of linear momentum equa- 
tion, the Cahn-Hilliard equation and the heat equation; the first two are the 
evolution equations for the two order parameters; the heat equation accounts 
for thermal effects like cooling or heat released due to a phase change. The 
model is complemented by a proper description of the free energy and con- 
sistent boundary and initial conditions; due to the presence of the gradients 
of the order parameters, two higher order boundary conditions appear to be 
necessary, and this can be also observed in the exploitation of the Clausius- 
Duhem inequality. 
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